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Abstract 

We present studies of the dc magnetization of thin platelike samples of the 
isotropic type II superconductor PbTl(10%), as a function of the angle be- 
tween the normal to the sample and the applied magnetic field H. We deter- 
mine the magnetization vector M by measuring the components both parallel 
and normal to H in a SQUID magnetometer, and we further decompose it 
in its reversible and irreversible contributions. The behavior of the reversible 
magnetization is well understood in terms of minimization of the free energy 
taking into account geometrical effects. In the mixed state at low fields, the 
dominant effect is the line energy gained by shortening the vortices, thus the 
flux lines are almost normal to the sample surface. Due to the geometrical 
constrain, the irreversible magnetization Mj rr remains locked to the sample 
normal over a wide range of fields and orientations, as already known. We 
show that in order to undestand the angle and field dependence of the modu- 
lus of Mj rr , which is a measure of the vortex pinning, and to correctly extract 
the field dependent critical current density, the knowledge of the modulus and 
orientation of the induction field B is required. 
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I. INTRODUCTION 



It is well known that sample shape effects have a strong influence in the magnetic response 
of superconductors. The discovery of the high temperature superconductors (HTSC) revived 
the interest in these effects, particularly for the case of platelike geometries |I]-[7j because 
those materials are frequently produced in the form of thin films or platelike single crystals. 

When we consider the shape effects on the dc magnetization of samples of type II su- 
perconductors whose dimensions are large as compared to the superconducting penetration 
length A, it is useful to divide the problem in two parts, namely the reversible magnetization 
and the irreversible magnetization that arises from vortex pinning. The first case appears as 
a rather simple textbook problem (at least for geometries that can be approximated by an 
ellipsoid) that can be solved by using the demagnetizing tensor formalism []8]|| . The second 
aspect is more complex, and have received considerable recent attention, both theoretically 
and experimentally ■ 

The macroscopic magnetic response due to flux pinning is well described by the criti- 
cal state model |10|,11|, according to which the current density flowing in any part of the 
sample is either the critical current density J c or zero. This current produces an spatially 
inhomogeneous magnetic induction B (r), that is detected in magnetometry studies as the 
non-equilibrium part of the magnetization. In general, J c is a function of the temperature 



T and the local B (r). The simplest critical state model, proposed by Bean flTT 



assumes 



that J c is independent of B, thus resulting in a uniform current density when the applied 
magnetic field H is varied in a typical isothermal measurement M(H) (hysteresis loop). 

Even within the simple Bean model, the calculation of B (r) in thin samples with H 
normal to the sample surface is rather complex, and has been the subject of extensive 
modeling and numerical analysis [|l|,0,|12j. The situation is further complicated when the 



realistic cases of field dependent and/or anisotropic J c (e.g., in HTSC) are considered |L3 
Several experimental techniques have also been used to measure the local B [|14],[3f 

Another aspect of the problem, that has been less explored, is the magnetic response 
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when H is tilted with respect to the normal to the surface of thin samples. As the non- 
equilibrium screening currents are strongly constrained by the geometry to flow parallel to 
the sample surface, there is a large angular range of applied fields in which the irreversible 
magnetization M irr points almost perpendicular to the surface |6]|7|], and consequently B 
is not parallel to H. This geometrical anisotropy results in some angular effects that are 
qualitatively similar to those observed in extremely anisotropic (quasi two-dimensional) 
superconductors |15|-|17|. The understanding of angular effects due to sample geometry are 



thus necessary for a proper interpretation of those studies. 

The above considerations indicate that it makes sense to explore the angular behavior 
of thin samples of isotropic superconductors. This is the approach followed, for instance, by 
Hellman et al. ||, and more recently by Zhukov et al. [[5J. Those studies are mostly focused 
on the angular dependence of the irreversible magnetization, and little attention is paid to 
the modulus and angle of the magnetic induction B, that represents the actual density and 
direction of the vortices. 

In this work we show that in order to understand the angle and field dependence of the 
irreversible magnetization of thin samples, and to correctly extract the field dependent J c 
using the critical state model, a complete knowledge of the vector B is required. We present 
studies of the dc magnetization of a platelike isotropic (conventional) superconductor as 
a function of angle with respect to H. We measure the total magnetization vector M and 
decompose it in its reversible and irreversible parts. We demonstrate that they have different 
directions, and analyze the angular and field dependence of both of them. We show that 
the geometrical effects on the reversible response determine the vector B, which is the key 
variable in the problem. 

II. EXPERIMENTAL DETAILS AND PROCEDURES 

In the present study we used samples of the isotropic superconducting alloy P&o.9^o.i- 
Two samples, a circle of 3.3mm in diameter and 43/im thick and a square of 2.4 x 2.4mm 2 
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and 39/im thick were cut from the same foil. The thicknesses were determined from the 
samples' weight and known density (~ llAg/cm 3 ), and later confirmed from the magnetic 
measurements. Samples dimensions were chosen so that the aspect ratios are similar to that 
of the Y Ba2Cu^,Of-i crystals that we use. The critical temperature is T c ~ 7 .3K. As results 
from both samples are very similar, only data from the square sample will be presented here. 

Dc magnetization measurements were made in a Quantum Design magnetometer 
equipped with two SQUID detectors, each one coupled to a set of pick up coils. This 
allows us to record two perpendicular components of the magnetic moment, namely the 
longitudinal component mi (parallel to H), and the transverse component mt (normal to 
H), and thus to determine the total magnetic moment vector m. Samples can be rotated 



in situ around an axis perpendicular to H using a home-made sample rotator device p8 
The angular resolution, given by the minimum step, is ~ 0.5°, but the overall error in the 
absolute value of the angle is about ~ 2°. The square sample was mounted with one pair of 
its sides parallel to the rotation axis. 

The initial alignment was performed as follows [|Tj|. Sample was first cooled in zero field 
below T c and then a small field was applied. Sample was then rotated using the home-made 
system until the angle formed by its normal and H was about 45°. Subsequently, the sample 
was rotated around the vertical axis (parallel to H) and the transverse component m t was 
measured as a function of that orientation. Once the maximum in m t was found, indicating 



that m lies in the plane formed by H and the axis of the transverse coils ||19|| , the angular 
position around the vertical axis was maintained fixed throughout the rest of the experiment. 
(The capability to rotate around the vertical axis is provided by the commercial system and 
is not affected by our perpendicular rotation device). 

It is known that the measurement of m t posses some difficulties which originate in the 
appearance of an spurious signal due to the longitudinal component m h that is detected 
by the transverse pick-up coils. This occurs when the sample is slightly off-centered with 
respect to the vertical axis of the coils [fill, which is frequently the case. We have completely 



and satisfactorily solved this problem. The solution includes the external processing of the 



original SQUID output signal using a software developed ad-hoc. All the details related 
with the hardware and software of the sample rotation system are presented elsewhere |L8| . 

We performed measurements of both components of the dc magnetization, Mi = mi/V 
and M t = m t /V (where V is the sample volume), as a function of applied field, at fixed 
temperature and for fixed values of Oh, the angle formed by H and the normal to the 
sample surface, n (see sketch in Figure |l]). We repeated these isothermal hysteresis loops for 
different Oh- After each loop was finished, the sample was rotated to the new Oh, warmed 
up above T c and cooled down in zero field (ZFC) before starting the new run. In this way, 
the initial Meissner response was recorded for each angle. 

III. RESULTS AND DISCUSSION 
A. Meissner response as a function of angle 

Figure [l] shows a typical pair of curves of M t and Mi versus H for the square sample, 
recorded at T = 4.5K using the procedure described in the previous section. In this particu- 
lar case Oh — 60°. Both M t and Mi exhibit an initial linear dependence on H characteristic 
of the Meissner response, followed by an hysteretic behavior indicative of vortex pinning 
when H is further increased and then decreased. 

We first analyze the Meissner response using the standard demagnetizing tensor formal- 
ism |§J|. Let's consider a homogeneous magnetic material immersed in a uniform field H. 
If the shape of the sample can be approximated by an ellipsoid, then the internal field Hj, 
the magnetization M and the magnetic induction B = H; + 47rM are all uniform, and the 
internal field is Hj = H+47T77M, where V is the demagnetizing tensor. We thus have 

B = H + 47r(T — 77) M (1) 

(Note that Eq. [I] is valid only if M is uniform. It cannot be used, for instance, to describe 
the critical state of a superconductor). 
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We apply Eq. [I] to the Meissner state, where B = 0. We decompose H in the two 
directions n and p, where the unit vector p is defined by the intersection of the sample 
surface and the plane formed by n and H (see sketch in Figure [I]). Then V becomes diagonal 
|nj@| and we have H\\ = HsmB H = -4tt(1 -v\\)M\\ and H± = HcosB H = -4n(l-v±)M±. 
For very thin samples u\\ ~ v and u± « 1 — 1v with !/<l, thus 

47r mi = _ (— cos 2 Q H H — sin 2 9^ VH 

\2i> 1 — 1/ / 

47rm t = - f ^ — \ sin ® H cos® H VH (2) 

\2i/ 1 — z/y 



where we have written eqs. [| in terms of the measured quantities m; and m t . A simple 
way to test this dependence is to make a ZFC, then apply a small field (typically ~ lOOe) 
and rotate the sample at fixed H . One of such measurements is shown in Figure [| where 
the solid lines are fits to Eqs. §. The small discrepancies between the data and the fits 
are due to the fact that the ZFC procedure is not perfect and thus there is always a small 
remanent field. This field induces a remanent magnetic moment in the sample that adds 
to the Meissner signal in Eqs. ^|. The periodicity of this remanent moment is 360°, instead 
of the 180° periodicity of the Meissner signal. When that effect is taken into account, we 
obtain excellent fits to the data in fig. [| as described elsewhere Jl8 . 



A more exact method is to determine the longitudinal and transverse Meissner slopes 
m 'i = ~~m an d m t = int fr° m t ne initial part of the magnetization loops M(H). This 
procedure eliminates the uncertainty in the value of H due to the remanent field, and also 
that arising from remanent magnetic moments, as they are field independent. In the inset 
of figure we have plotted m[ and m' t as a function of Qh for all the angles measured. 
We obtain excellent agreement with eqs. ^| with the parameters V = 2.29 x 10 _4 cm 3 and 
v = 0.018. The volume so obtained coincides within 1% with the volume calculated by 
weighting, and it is the value that we have used to calculate the curves of figure ^ and all 
the magnetization data shown from now on. The demagnetizing factor coincides very well 
with the expected value v ~ t/L = 0.016. 
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It is also useful to calculate the modulus and orientation of M in the Meissner state from 
Eqs. |2|. If we call a the angle formed by M and H, then tana = (m^/mj) is 

'1 - 2v) tan9# 



tana 



1 + 2v tan 6^ 



(3) 



while the modulus is 



- 4vrM 



sin 9^ 



cos 6 



2 -|V2 



U 



H 



(4) 



[l-v) 2 (2z/) 2 

We found that Eqs. ||] and [3] allow us to determine the angle 0# from the initial slopes 
of our M (H) data with accuracy better than 1°. As this error is lower than that associated 
with our mechanical devise (~ 2°), we decided to use in all cases the value of 0# obtained 
from the Meissner slopes, and take the reading of the rotator as a double-check. 



B. The reversible response in the mixed state 

In this section we will show that the reversible response of the sample behaves as expected 
in an isotropic superconductor. Assuming that the Bean model is valid, both the longitudinal 
and transverse components of the equilibrium or reversible magnetization M e(J in the mixed 
state, M eq j and M eg>t , can be determined by simply taking the mean value between the 
two branches of the hysteresis loops, as illustrated in figure [1} We see that for fields above 
~ 500Oe the transverse component M eqjt ~ 0. This is true for all 0#, i.e., at high fields 
M eg is parallel to H regardless of the sample orientation. 

In figure |3] we show the modulus M eq as a function of H for several field orientations. 
Also shown are the corresponding Meissner responses, which are angle dependent as expected 
from Eq. f|. On the other hand, in the mixed state at high fields (above ~ 500Oe) all the 
curves coincide; and in particular the condition M eq « that determines the upper critical 
field H C 2 is independent of 6#. 

The behavior of M e(? in the mixed state at high fields is a consequence of the isotropy of 
the material. Indeed, the free energy F of the vortex system in isotropic superconductors 
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is only a function of B, and the equilibrium state is obtained by minimization of the Gibbs 
free energy [see Eq. 8.60 in fllEfl l 



at high fields 4irM eq <C H, then B«H, and the last term in Eq. |5| is negligible. We thus 
recover the solution of an infinite superconductor, where B, H and M eg are parallel and 
M eq (H) is given by the usual expressions f| regardless of the sample geometry. 

The penetration of the first vortex becomes energetically favorable when the modulus 
of the internal field, Hi, equals the lower critical field H c \. Taking into account that in the 
Meissner state H = (1 — 17) Hi, we can define an "apparent" lower critical field H* x (Oh) 
that, in spite of the material isotropy, is angle dependent due to geometrical effects: 



if*, = H c i 



/"v 2 o 21 -1/2 

sin <c>h _j_ cos Hfl 



(6) 



The experimental determination of H* x (0#) is complicated by the fact that for fields 
only slightly above the Meissner regime, the critical state is not fully developed in the 
increasing branch of the loop (see figure Q), and the mean magnetization does not give us a 
good estimate of M eq . To solve this difficulty we adjusted the M eq (if) data at intermediate 
fields by a dependence M eq oc In if. Although this logarithmic dependence is only expected 
in high k superconductors |T6|,|8||, while in the present case k ~ 2.2, the experimental fact is 



that the fit works pretty well. The upper limit of the field range of validity of the fit is about 
H ~ 500Oe. The lower limit depends on Oh, varying from H ~ 150Oe for 0# ~ 0° to 
H ~ 250Oe for 0# approaching 90°. We then identified the field at which the extrapolation 
of the fits intersect the Meissner slope with (0#), as shown in figure [| 

Evidence that the extrapolations provide good estimates of M eq come from the observa- 
tion that the intersection with the Meissner slopes occur at approximately the same mag- 
netization value AirM eq ~ 290Oe for all angles. This is what we would expect, as from Eqs. 
[| and [] results 47rM eg (if* 1 ) = H c ±, independent of 0#. 
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C. The vortex orientation 



We can now use the M eg data to calculate the equilibrium induction field B eg according 
to Eq. |l|. In figure [| we plot the angle 9 b between B e? and n, as a function of H for several 
orientations 9#. At high fields Ob ~ Qh, as already discussed, but as H decreases 0# 
approaches to zero for all orientations. 

This behavior is also understood from thermodynamic considerations. From Eq. ^ we 
can obtain the equilibrium condition by minimization with respect to B eg , i.e. d ^ G = 



do 



eq, 



0. Combining with Eq. [I] these conditions result 
dF B eq> ± {B eq ,i_ — H±) 



OF B eq \\ (B eq> \\ - H\ 







dB eqA 4vr 4?r (1 - v) 



(7) 



For !/<l, the first expression implies that B eqt ± ~ H± and the second condition reduces 
to d g F - ~ For fields slightly above H* 1: vortex density is very low (B eq H c i), then 
interactions are negligible and F ~ (B eq /<&o) Ei, where E\ is the vortex line energy [i~6||. In 



this limit B eq \\/H\\ ~ B/H c i, thus B eq \\ Hn, i.e., the vortex direction is very close to n. 
At intermediate fields F w ^ + ®° B *? 2 In (H c2 /B eo ) and we obtain 

07T 2(47rA) h 

B,„^m-^^J H ^ 



The physical interpretation of these results is straightforward. The normal component 
remains almost unchanged due to flux conservation (a result that of course becomes exact 
in the limit of an infinite slab). At low fields the parallel component is reduced because 
the system gains energy by shortening the vortex length. As field increases interactions also 
favor the reduction of B eq n, but eventually the energy cost of having B eg non-parallel to 
H becomes too large. Thus, vortex orientation shifts from n to H as field increases. The 
modulus B eq evolves from B eq ~ H± at low fields to B eq ~ H at high fields. 
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D. The irreversible response 



According to the Bean model [11], the irreversible magnetization is related to J c through 
a geometry dependent proportionality constant. We define the irreversible magnetization 
vector as Mj rr = | AM, where the components of AM are the widths of the longitudinal and 
transverse hysteresis loops, AM; and AM t respectively. In thin samples the non-equilibrium 
currents are constrained to the sample plane, thus they should generate an irreversible 
magnetization Mj rr that is almost locked to the sample normal J/J in a wide angular regime 
< Qh < @o where tan (6 C ) ~ L/t. When Q H ~ 9 C , the direction of M irr is expected 
to rotate ~ 90° in a narrow angular range and to locate in the plane of the plate. This 
behavior has been recently demonstrated by Zhukov et al. @ in detailed studies of the 
modulus and orientation of the irreversible magnetization of both isotropic and anisotropic 
thin superconductors. 

We have confirmed that in our samples, for 0# < 70° and for all applied fields, M irr 
points in the direction of n within our 1° resolution. For angles G# > 70° the situation is 
more complex. Mj rr progressively tilts away from n as B# increases, but its direction is 
field dependent. This behavior cannot be described by the expressions proposed by Zhukov 
et al. 0, where the field dependence of J c is not considered. 

In figure |5](a) we show the modulus M irr as a function of H for different angles 6#. 
For clarity we only show data up to 0# slightly above 70°, thus all the curves in the figure 
correspond to M irr || n or very close to it. We first note that M irr goes to zero for all Qh 
at H = H C 2, as expected for an isotropic type-II superconductor. Another evident result is 
that, although our superconductor is isotropic, M irr is angle dependent, even for Q H < 70°. 
We will now show that this dependence originates in sample geometry effects. 

Clearly, to understand these results we must go beyond the Bean model and consider 
the field dependence of J c . There are two basic points to be taken into account. The 
first one is that J c (and consequently M{ rr ) is a function of B rather than H. The second 
one is that the non-equilibrium currents determine the direction of the variation of B (r) 
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(through J oc V x B) but not the vortex orientation itself. Indeed, the field distribution in 
the critical state of thin samples, either for 0# = or for 0# 7^ 0, has been calculated 
taking into account the contribution of the irreversible currents plus the applied field, but 
not the reversible magnetization []]]. 

We will consider that the field B (r) arises from the superposition of the spatially in- 
homogeneous contribution produced by the non-equilibrium currents and the homogeneous 
reversible contribution calculated in the previous section. Provided that the irreversible 
contribution to B (r) is not too large as compared to B e(? , as indeed occurs in our samples, 
then the average vortex orientation is still approximately given by B eg . This also allows 
us to make the approximation that M irr is simply a function of B eq . (This is equivalent to 
the usual approach of considering that, although J c is field dependent, it is almost constant 
over the range of field variation within the sample for a given H, and thus the Bean model 
can be applied to each field). At low fields we can further approximate B eq ~ H±, thus we 
expect M irr to be a function of H_]_. 

To check the above modeling, in figure |^(b) we re-plot the M irr data of fig. [5](a), as a 
function of H±. The good scaling of the curves in the low field range indicates that the 
angular dependence of M irr in this isotropic superconductor can be understood in terms of 
the angular-dependent equilibrium vortex orientation. 

In figure^ we show M irr as a function of for several H. At very low fields (H < lOOe), 
M irr is almost constant over a wide angular range, starts to decrease at around 0# ~ 70° 
and reaches its minimum value at G# = 90°. This is, at least qualitatively, the expected 
behavior [[?[] in thin samples when the field dependence of J c is not taken into account. 
Indeed, the dotted line is the dependence predicted by Zhukov et al , where we have used 
the aspect ratio v ~ 0.018 of our sample determined previously, and thus there are no free 
parameters. We note, however, that a better fit could be obtained using a larger v. 

At higher fields we must take into account the vortex orientation as described above. On 
the other hand, for 0# < 70° we do not need to consider the reduction of M^rr due to the 
geometrical effects described in [fjj. We thus expect an angular dependence 
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M irr (H,Q H ) = M irr (H ± ,Q H = 0) = M° r (# cos Q H ) (8) 

where Mf rr is the curve measured for Q H = 0. The validity of this scaling is manifested 
in the figure []. An excellent agreement between the experimental data and eq. |] is observed 
for Qh < 70° in the field range < H < 240Oe, where B is approximately aligned with 
the normal to the sample. For 0# > 70° the decrease of Mj rr due to sample geometry 
takes over. For fields above 240Oe, the fit fails because the B eq ~ H± approximation is not 
satisfied. We note that this model satisfactorily accounts for the angular dependence of M irr 
that Zhukov et al. observed at high fields in their isotropic samples, which could not be 
explained within their picture. 



IV. CONCLUSION 

We have shown that the variable that governs the angle dependence of M irr in isotropic 
thin samples is the magnetic induction B, which describes the density and orientation of the 
vortices. The behavior of M irr (H, 0#) at low fields and over a wide angular range can be 
understood by simultaneously taking into account the modulus and orientation of the equi- 
librium magnetization vector B e(? , and the field dependence of Mj rr when H is normal to the 
sample. To perform the analysis reported here, it was necessary to simultaneously measure 
both components of the magnetization in an isotropic sample of well defined geometry, and 
to be able to decompose it in the reversible and irreversible parts. This last step was only 
possible because we used a sample with very low J c . 



It is a well known experimental fact p0[ that HTSCs exhibit the same scaling of Mj, 



with H± as described by Eq. ||. That behavior is consistent with the anisotropic scaling 



expected for these materials in the limit of very large anisotropy |15|JT6|| , and is usually taken 
as evidence for quasi two dimensional behavior. However, we have observed that the same 
scaling occurs in our isotropic sample, thus indicating that geometrical effects should be 
carefully taken into account in those studies. 



12 



REFERENCES 

[1] Clem J R and Sanchez A 1994 Phys. Rev. B. 50 9355 

[2] Brandt E H 1994 Phys. Rev. B. 49 9024; 1994 Phys. Rev. B. 50 4034; 1995 Phys. Rev. 
B. 52 15442 

[3] Schuster T, Kuhn H, Koblischka M R, Theuss H, Kronmiiller H, Leghissa M, Kraus M 
and Saemann-Ischenko G 1993 Phys. Rev. B. 47 373; Schuster Th, Kuhn H, Brandt 
E H, Indenbom M V, Klaser M, Miiller-Vogt G, Habermeier H-U, Kronmiiller H and 
Forkl A 1995 Phys. Rev. B. 52 10375; Schuster T, Leghissa M, Koblischka M R, Kuhn 
H, Kraus M, Kronmiiller H and Saemann-Ischenko G 1992 Physica C 203 203 

[4] Frankel D J 1979 J. Appl. Phys. 50 5402 

[5] Kunchur M N and Poon S J 1991 Phys. Rev. B. 43 2916 

[6] Hellman F, Gyorgy E M and Dynes R C 1992 Phys. Rev. Lett. 68 867 

[7] Zhukov A A, Perkins G K, Bugoslavsky Yu V and Caplin A D 1997 Phys. Rev. B. 56 
2809 

[8] Parks R D (ed.) 1969 Superconductivity Marcel Dekker Inc., New York 

[9] Landau L D and Lifshitz E M 1960 Electrodynamics of Continuous Media Pergamon 
Press Ltd., Oxford 

[10] Campbell A M and Evetts J E 1972 Adv. Phys. 21 199 

[11] Bean C P 1962 Phys. Rev. Lett. 8 250 

[12] Daeumling M and Larbalestier D C 1989 Phys. Rev. B. 40 9350 
[13] Conner L W and Malozemoff A P 1991 Phys. Rev. B. 43 402 

[14] Tamegai T, Krusin-Elbaum L, Santhanam P, Brady M J, Feild C and Holtzberg F 1992 
Phys. Rev. B. 45 2589 

13 



[15] Blatter G, Geshkenbein V B and Larkin A I 1992 Phys. Rev. B. 6 875 

[16] Blatter G, Feigel'man M V, Geshkenbein V B, Larkin A I and Vinokur V M 1994 Rev. 
Mod. Phys. 4 1125 

[17] Buzdin A I and Simonov A Yu 1991 Physica C 175 143 

[18] Casa D M et al. to be published 

[19] Quantum Design operating manual and application notes 

[20] Roas B, Schultz L and Saemann-Ischenko G, 1990 Phys. Rev. Lett. 64 479; Schimtt P, 
Kummeth P, Schulz L and Saemann-Ischenko G 1991 Phys. Rev. Lett. 67 267 



14 



FIGURES 

FIG. 1. (a) Open symbols: Transverse component of the magnetization as a function of applied 
field for the sample normal at an angle of ~ 60°. The arrows indicate the direction of the field 
sweep. Full symbols: Mean value of the magnetization, (b) Idem for the longitudinal component. 
The sketch shows the components and angles used in the text. 

FIG. 2. Angular dependence of the longitudinal and transverse components of the magnetiza- 
tion in the Meissner state, as measured by rotating the sample at fixed applied field (the origin of 
angles is arbitrary). The solid lines are fits to Eq. S Inset: Meissner slopes as a function of angle. 
The solid lines are fits to Eq. ||, see text. 

FIG. 3. Modulus of the equilibrium magnetization in the Meissner and mixed states as a 
function of applied field, for several sample orientations. The dotted lines are fits M eq oc InH. 
The inset shows a blow up of the same data in semilogarithmic scale. The intersections of the 
extrapolations of the logarithmic fits with the Meissner response (shown by crosses) provide an 
estimate of H* t (Oh) 

FIG. 4. Angle between the equilibrium induction field B and the normal to the sample, as a 
function of applied field, for three sample orientations Oh- The symbols are experimental data; 
the dotted lines are the extrapolations of the logarithmic fits shown in figure 3. 

FIG. 5. (a): Irreversible magnetization as a function of applied field for several sample orien- 
tations, (b): Same data, as a function of the normal component of the applied field. 

FIG. 6. Irreversible magnetization as a function of Oh for several applied fields. From top to 
bottom, H = 0, 20, 30, 40, 50, 60, 70, 80, 100, 120, 160, 200, 240Oe. The solid lines are fits to Eq. |. 
The dotted line is a fit to the model of Ref. |7j using an aspect ratio u = 0.018 
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